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Abstract 

Correlated Basis Function perturbation theory is used to evaluate the zero 
temperature response S{q,uj) of ^He-^He mixtures for inelastic neutron scat- 
tering, at momentum transfers q ranging from 1.1 to 1.7 A^^. We adopt a 
Jastrow correlated ground state and a basis of correlated particle-hole and 
phonon states. We insert correlated one particle-one hole and one-phonon 
states to compute the second order response. The decay of the one-phonon 
states into two-phonon states is accounted for in boson-boson approximation. 
The full response is splitted into three partial components Sai3{q,io), each of 
them showing a particle-hole bump and a one phonon, delta shaped peak, 
which stays separated from the multiphonon background. The cross term 
S34:{q,uj) results to be of comparable importance to S'33(g,i-i;) in the particle- 
hole sector and to 5*44 (g, i^) in the phonon one. Once the one-phonon peak has 
been convoluted with the experimental broadening, the computed scattering 
function is in semiquantitative agreement with recent experimental measure- 
ments. 

67.55.Lf, 67.60.-g 

Typeset using REVTgX 

1 



1. INTRODUCTION 

Isotopic atomic Helium mixtures are an intriguing case for many-body physicists. There 
exists a large body of experimental data, concerning mostly static properties (for instance, 



quantities, as the zero concentration {x^ = 0) ^He effective mass (^3), have been also 
measured^. Recently, inelastic neutron scattering experiments have been carried out both 
at low, or intermediate, and high momentum transfersifl. In both regions the measured 
response presents two generally distinguishable structures, to be ascribed to boson like col- 
lective excitations (phonons and rotons) and to Fermi particle-hole ones. However, this 
apparently simple picture hides a large interplay between the components of the mixtures, 
each of them probably contributing on comparable foot to both the branches of the re- 
sponse. The reason for this lies in the large correlation effects, which are present in the 
system because of the strong interatomic potential and of the large density. This are also 
the motivations why truly microscopic and ab initio studies of Helium mixtures are difficult, 
and, in the case of the response, practically absent in literature. 

Qualitative studies of the response have been done in ref.i (using a matrix dispersion- 
relation representation) and in ref.i with a correlated RPA approach (very similar, in spirit, 
to the phenomenological Polarization Potential method used in ref.0. Here we will employ 
the Correlated Basis Function (CBF) perturbation theory, to embody the above correlation 
effects directly into the basis functions. CBF has shown to be a powerful tool to succesfuUy 
study Helium at zero temperature: the energetics of both pure ^He and '^He are well described 
by sophisticated correlated ground state wave functions, containing explicit two-, three- 
body, back-flow and spin correlation^lHEi; properties of one ^He impurity in ^He, such as 
chemical potential and effective mass are also quantitatively reproduced by such correlated 
wave functions^. In particular, by using CBF based perturbation theory, with the insertion 
of up to two correlated independent phonon intermediate states, the impurity effective mass 
mg turns out to be 2.2m3, to be compared with the experimentally measured 2.3m3 value. 

The behavior of the ^He effective mass with the concentration in dilute mixtures has 



the chemical potentials and the maximum 




Excitation spectra and related 
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been recently object of some debate. Specific heat measurement slll'lll at finite X3 do not show 
appreciable deviations from its = value. In ref.l the Authors have to postulate a much 
larger value (mg ~ 2.95m3 at 0:3 = 0.05) in order to reproduce the position of the particle- 
hole response with a Lindhard like function and using a simple Landau-Pomeranchuk (LP) 
quasiparticle spectrum0, 

This contradiction does not appear if one modifies the LP spectrum ( LP modified, or 
LPM) as: 

e,{LPM)=eo + ' (2) 

There are both experimental and theoretical! indications of a deviation from the simple 
LP form. 

In a CBF based approach, we assume to have an homogeneous mixture of A'"3 ^He atoms 
and A'4 ^He atoms in a volume f2, with partial densities Pa=?,,i = N^/^, total density 
p = Ps + Pi and concentrations Xa = Pa/ P- We will keep constant densities, while letting 
Na and the volume going to infinity. The nonrelativistic Hamiltonian of the mixture is 

^ = - E E^v? + ^ E E n-.), (3) 

a=3,4 i=l ^'"-a ^ a,/3=3,4 i^j 

where the interaction is the same for all the different pairs of the mixture. 
A realistic, correlated, variational ground state wave function \l'o is obtained by the 
Jastrow-Feenberg ansatzEl 

^0 = FjFtFbfMNs) , (4) 

where (poi^s) is the ground state Fermi gas wave function for the ^He component and Fj, Ft 
and Fbf are N-body correlation operators including explicit two-, three-body and back- 
flow dynamical correlations respectively. We will limit our analysis to the case of two-body, 
state independent (or Jastrow) correlations only {Ft = Fbf = Fj results to be 



N3 N4, Ns N4 

FAN„m) = n f'''\nsn) n f^'''\n.u) nn/^'-'H-^aj- (5) 

«3<j3 «4<j4 «3 ii 

where /(°'^)(r) are two-body correlation functions determined by minimizing the variational 
ground state energy. 

It is possible to generate a correlated basis through the operator (j^), to be used in a CBF 
perturbation theory (CBFPT). This theory has been succesfuUy adopted for computing the 
inclusive response of nuclear matter and heavy nuclei to electron and hadron scattering&0, 
and has shown to be able to provide a semiquantitative agreement with experimental neutron 
inelastic scattering (nIS) data in pure, liquid atomic ^Heil. 

In this paper, we will apply CBFPT to compute the nIS response of the mixture, by 
considering as intermediate states the normalized, correlated ^He n-phonon states (nPH) 
|ki,..,k„), and ^He n-particle, m-hole states (np-mh) |pi, .., p„, hi, .., h^, )• 

The nPH states are given by 

Ij^ P4(ki)..p4(kn)|^o) 

(vl/o|pl(k„)..pi(ki)p4(ki)..p4(k.)|v^o)^/^' ^' 
where P4(k) is the ^He density fiuctuaction operator 

P4(k) = E e^'-""' • (7) 

Correlated np-mh states are obtained in a similar way, by applying the correlation op- 
erator to the Fermi gas excited states ^np-mhiNs), 

Ipi, .., p„, hi, .., h^) = — . (8) 

\'^np-mh\-f j-f j\'^np-mh) ' 

We will consider 1-phonon (IPH) and Ip-lh intermediate correlated states, which we 
will term as One Intermediate Excitation (OIE) states. The response computed at the OIE 
level will be called variational. In addition, we will also consider the possible decay of IPH 
states into 2PH ones, which is essential in giving a physically meaningful ^He excitation 
spectrum and provides a quenching of the one-phonon peak. This term will be computed 
in a boson-boson approximation, i.e. neglecting the ^He antisymmetry. Such an approach 
may be justified on the basis of the low ^He concentration. 



Ip-lh states may also be coupled to IPH and 2PH. Such a coupling may be taken 
into account by a corresponding self-energy insertion. Its analogous in the problem of the 
single ^He atom in ^He is responsible for the impurity large effective mass. To estimate the 
importance of this effect we will use the on-shell part of the impurity self-energy, again 
relying on the small value of X3. 

The plan of the paper is as follows. In section II we will briefly outline the CBFPT 
for the response of the mixture and the variational calculation will be described in some 
details. Section III is devoted to the description of the calculation of the coupling with the 
2PII states and of the decays into IPH and 2PII states. Section IV contains results for the 
response and the comparison with the experimental scattering functions. Moreover, the ^He 
and ^He excitation spectra are presented and discussed. Conclusions are drawn in section 
V. 

2. CBFPT FOR THE RESPONSE 

The Dynamical Structure Function (DSF) S{q,u!) of a ^He-^He mixture at T = is 
given by the imaginary part of the polarization propagator D{q,uj) 

S{q,u;) = -QD{q,u;), (9) 



where 



D{q,uj) = l(^o|pt(q)^ -p(q)|^o), (10) 



and 



P(q) = P3(q)+P4(q), (11) 

Pa(q)= E e^'*-'"' (12) 

i=l,Na 

and N = + A^4. In eg. ([TI]|) , is the exact ground state of H with eigenvalue Eq. 
The total DSF may be expressed in terms of partial a/5 DSF, Sai3{q,uj), as 

S{q,uj)= J2 XapSapiq.uj), (13) 

a,/3=3,4 



with Xai3 = {xaXpY^"^ and 

= \^D^,(,,.) = i3-^(*„|p.(q)^^-^^L_p,(q)|*„). (14) 

The experimentally measured nIS double differential cross section for the mixture directly 
provides access to the the total scattering function S{q,u), wich is in turn related to the 
partial DSFs' by the relation 



The elementary cross sections, as given by Sears" in units of barns, are (T4 = 1.34, cxa = 
4.42, a\ = 1.19 and as4, = 4.70. The inchoerent DSF 5*33(5, tu) also appears in the expression 
above, with the corresponding cross section a^. Because of the absence of spin correlations, 
it describes ^He spin fluctuactions via the operator 

P3(q)= E e^'^-^'i^^ (16) 

i=l,N3 

where Ij is the spin of ^He i-nucleus. 

We will focus, in the remainder, mainly on the calculation of Saf3- To derive a perturba- 
tive expansion it is convenient to split H into an unperturbed piece Hq and an interaction 
term Hi, as follows 

{m\Ho\n) = 6nm{m\H\m) = El^, (17) 

and 

{m\Hi\n) = (1 - 5nrn){m\H\n) = Hmn- (18) 

Here \m) are correlated basis states, eigenstates of Hq. In particular, |0) = |^o) is not 
an eigenstate of H and its difference from |\E'o) is treated perturbatively. The expansion is 
obtained by writing 

H-Eo = Ho-E^, + {Hi - AEo), (19) 
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where AEq is the correction to the variational ground state energy Eq, and by developing the 
propagator [H — Eq — ui — ir))"^ in powers of {Hi — AEq). A similar expansion is performed 
for the ground state |^o)- 

If the expansion is truncated at the zeroth order, the partial DSF are given by: 

Sapiq,u) = _^^(*o|pl(q)|n)(n|p;3(q)|*o>%-c^n), (20) 



with oJn^E-^- El 



As stated in the introduction, we will first consider only OIE insertions, i.e. correlated 
IPH and Ip-lh intermediate states, defined as: 

P4(k)|*o) 



(*o|pi(k)p4(k)|*o)V^ ' 



Fj\ 




-Ih) 




\FlFj\ 





(21) 



p,h) = ^J\J^p-^hi ^22) 



2.1 THE VARIATIONAL RESPONSES 

The variational response is given by the sum of two components, 

S^^{q,u;)^Sl$''{q,u;) + S'j'f'\q,u;), (23) 
where S]^^{q,uj) has a IPH intermediate state 

S'j;"{qM = ^=i=E(*o|pl(q)|k)(k|p^(q)|*o)%-c.,), (24) 

and S^p~^^{q,u!) has a Ip-lh intermediate state 

Sl'i'\q,u^) = ^=i=E(*o|pl(q)|p,h)(p,h|p;3(q)|^o)5(o;-ep + e.). (25) 

\JNaNp p^h 

cuk and Sp — Eh are the variational energies of the OIE states considered. 
ujk is given by 

_ 1 (k|if-Eg|k) _ h'e 
' N^ (k|k) 2m4Su{ky ^ ' 
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and corresponds to the well known Feynman spectrumB. In eg. (pB]) , 6*44 (fc) is the variational 
estimate of the 44 component of the Static Structure Function (SSF), Sapik), given by: 

1 (^o|pt(k)p^(k)|^o) 
^ap[li) - I /,T, |,T, \ • 



In a similar way, ex=p^h is obtained by 

(28) 



(x| 


\H - 




|x) 


(x| 


|x) 



where |x) is a particle or hole correlated state. We will discuss later the evaluation of e^- 
By using the definition of the SSF given in eq. (p7|) , ^Q,(q;k) = (\E'o|Pa(q)|k) is readily 
obtained as 

e„(q;k) = v/iV;^^£L4-q, (29) 

giving, for S'^J^(g,u;), 

-) = E - u,). (30) 

The one-phonon contribution to the variational a-f3 responses shows a delta-like be- 
havior, whose strenght is Z^^fjik) = Sai{k)Sfii{k) / Si^ik), and it is located at the Feynman 
phonon energy. We notice that (i) Zl^{k) = Su{k) and that (ii) the 33 and 44 variational 
DSF are positive (S'44(A;) being positive), whereas this may not be true for the 34 DSF. 

The expression for the particle-hole response S^^^^{q,uj) is more involved. A detailed 
description for a pure Fermi system (specifically, nuclear matter) can be found inlli and 
references therein. On the basis of that formalism, the extension to a boson-fermion mixture 
is straightforward. 

In CBF theory, the non diagonal matrix elements ^a{^] P? h) = (\I'o|pQ(q)|p, h) are com- 
puted by a cluster expansion in Mayer like diagrams, and by summing infinite classes of 
relevant terms. The C,a are explicitely given by: 

^a(q; p, h) = (5q-p+h^= ^ {hdd,a3{q) + (^asfl + KdMq)])^ (31) 

/D{p)D{h) 
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where 

hxy,a3{q) = P» J d^re^'^'^'lgxy^asir) - S^y^dd], (32) 

with (x, y) = {d, e) and gxy,a(3{'^) are partial radial distribution functions (RDF). In fact, the 
total a/5-RDF, (^^^(r), giving the probability of finding a a-type particle 1 at a distance ri2 
from a /3-type particle 2, 



iVa(iV/3-5a/3)^ 






i^or 




1 rfVi. 




i^or 



^a/3(ri2) = 7-- -, (33) 



is computed, in Fermi Hypernetted Chain (FHNC)li, using the correlated g.s. and it 
turns out to be written as: 

ga/dir 12) = 9dd,cp{r 12) + ^pz9de,a3{ri2) + 5 aZ9ed,Zp{r 12) + (5a3^/335'ee,33('^12)- (34) 

The partial RDF are classified, in FHNC theory, according to whether the external 
particle (1 or 2) is reached by a statistical correlation (i.e. if the particle is involved in an 
exchange loop, e-vertex), or by a dynamical correlation, {f^^'^^'Y ~ ^i only (rf-vertex). 

The definitions of the partial RDF, together with the full set of the related FHNC 
equations, may be found in@. 

Actually, eq. ( PUf ) sums all cluster diagrams factorizable in products of dressed, two-body 
diagrams. They do not contain only two-body cluster terms, but include, in turn, an infinite 
number of particles, as they are written in terms of the RDFs, rather than the bare two-body 
correlations. 

Three-body, non factorizable diagrams are also present in the cluster expansion of ^q, 
even if they do not appear in eq . (^91) . However, they have been inserted, following ref.lll. 
The function D{x = p, h) is 

D{x) = 1 - P3 / t^'re*""(^?dd,33(r) - l)L(A;^r), (35) 

where L{kFr) is the FHNC generalization of the exchange Slater function l{kFr) = 
?)ji{kpr) / {kprY , and kp is the ^He Fermi momentum {kp = Svr^ps). 

9 



Again, as D{x) turns out to be positive, 5'33_44''(g, a;) are positive, whereas Sll~^'^{q,u) 
may be not. 

Finally, the spin fluctuaction matrix element, ^3(q;p, h) = (\l'o|p3^(q)|p, h), is simply 
given by: 

^3(q; P' h) = Vp+h^= "^ (36) 

lD{p)D{h) 



3. CORRELATED ONE AND TWO PHONON INTERMEDIATE STATES 

In this section we will first study the effect on the phonon responses of the insertion of 
orthogonal, correlated 2PH states: 

|kik2)o=(l-|k)(k|)|kik2), (37) 

where the 2PH states of eq.(^ have been orthogonalized to the IPH ones by a Gram-Schmidt 
procedure. 

2PH states influence the partial polarization propagators Dai3{q,uj) via the direct cou- 
pling to the ground state and via the decay of IPH states into 2PH. The coupling to the g.s 
goes through the matrix element of the ^He fluctuaction operator, 

e3(q;ki,k2) = (v&o|pkq)|ki,k2)o, (38) 

( notice that ^4(q; ki, k2) vanishes because of the Schmidt orthogonalization of the 2PH 
states), whereas the decay is driven by the non diagonal matrix element of the hamiltonian 

a(k;ki,k2) = {k\H,\k,M)o. (39) 

These CBF matrix elements have been computed in a boson-boson approximation (treat- 
ing the ^He as a mass-3 boson) and by adopting the Convolution Approximation (CA) for 
the three-body distribution functions. 

Their explicit expressions are: 

e3(q; k„ k2) = f"^^^^^"^^^^ f ^33(g) - , (40) 

/^44(A;i) ^44(^2) V 
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and 



fk-k k)- / k- ki5'44(fc2) + k ■ k25'44(fci) - fc^5'44(fci)5'44(A;2) 



[x^ k'^S34{k)S34{ki)S34{k2) 



^3 5*44 (/c) ^5*44 (/c) 5*44 (/Ci) 5*44 (/C2)^ 

It is convenient, at this point, to introduce the correlated self-energy 

= 1 J: If-'^-'^-'r . (42) 

and the function x(q; /c, cu) given by: 

x(q;A;,u;) = J-^ a(k;ki,k2) — -CUd-MM)- (43) 

If we define the dressed phonon propagator G'^{k, lo) as 

G'{k, = -— ^- -, (44) 

ook — ^1 [k, uj) — uj — irj 

then, the phonon contributions to the polarization propagators can be rearranged as: 

D^.rM = ^El(*o|pi(q)|k)rG'^(A;,a;), (45) 
JV4 k 

D^^^M = ^^E(*o|pi(q)|k)G'^(^,a;)[(k|p3(q)|*o) + ^x(q;^,^)], (46) 
and 

D'^i'{<l.^) = ^E(^o|4(q)|k)G'^(^,(^)[(k|p3(q)|*o) + -^x(q;^,^)]. (47) 

The DSF are then obtained by taking the imaginary parts of D^fj ■ Terms quadratic in 
^3(q; ki, k2) have not been considered. 

The relevant changes introduced by the insertion of the 2PH states in the phonon re- 
sponses are: 

1. the strengths of the delta-like IPH peaks Zap are generally quenched respect to Z^^. 
In the 44 case, we have 
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Zu(k)^ZUk)(l + ^^^£^y . (48) 



Analogous corrections occur for Z34{k) and Z^^lk), which are also affected by those 
parts of the polarization propagators containing ^3(q; ki, k2); 

2. the IPH peaks are shifted by the real part of the on-shell self energy, since the ^He 
spectrum is modified as 

uJk uj^""^ = oJk + 3?Si(A;, u'i^^)- (49) 

3. a multiphonon tail appears at large w-values, beyond the position of the IPH peak, 
at the momentum transfers here considered. 

4. CBF RESPONSES 

In the class of the Jastrow correlated wave functions, the best variational choice is pro- 
vided by the solution of the Euler equations 

The resulting equations have been derived, within the FHNC framework, and solved 
for the ^He impurity problemBU, for the boson-boson mixtureil and, lately, for the real 
fermion-boson casei. 

Another, often used approach consists in parametrizing the correlation functions and in 
minimizing the ground state energy with respect to the parameters. This is the choice we 
have adopted here. Besides that, some of the results we will present have been obtained 
within the Average Correlation Approximation (ACA). In ACA, the correlation functions 
are the same for all the types of pairs (/'■^'^'' = Z*-^''^'' = Z*-^''^'') and the differences in the 
distribution functions (or in the static structure functions) are due only to the different 
isotope densities and statistics. We will also show that going beyond the ACA does not 
affect our results. 
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We have used three types of correlation functions: the time honored, short ranged McMil- 
lan form (SR) and two long ranged functionsa (LR and LRl). 
The McMillan correlation, in ACA, is given by: 



fsR{r) = exp 



%a\^ 1 



(51) 



where b = 1.18 and a = 2.556A. The SR correlation function gives a good description 
of the short range behavior of the pair wave function but fails to reproduce long range 
properties. For istance, it does not ensure the linear behavior of the ^He SSF at /c (the 
phonon dispersion). Such a dispersion reflects in a long range behavior of the correlation of 
the type f{r^ oo) — 1 oc — r^^. To this aim, we have also used a modified form, having the 
correct long range structure (LR), given by 

fLR{r) = fsR{r) A + Bexpl^-^^^^^^^ . (52) 

The parameters of /Li?(^), giving the variational mimimum of the ^He energy at the 
equilibrium density po = 0.02185A-3, are h = 1.18, A = 0.85,5 = I- A,D = 3.8A and 
r = 0.043A~^ (see Ref.ii for more details about the energetics of the mixture). The B and 
r parameters are related to the experimental pure ^He sound velocity c and to the low-Zc 
behavior of its SSF by the relations: 

- = ^;S^'^Hk^^) = ^. (53) 

In order to check the accuracy of ACA, we have also used a LR correlation (LRl), 
formally identical to fiR, but with parameters depending on the type of the correlated pair. 
The 44 correlation function is the same as above, whereas the parameters of the 43 and 33 
ones have been obtained by minimizing the energy of the pure ^He with one and two ^He 
impurities, respectively@. 

Key ingredients in the CBF theory of the response in Helium mixtures are the radial 
distribution functions Qapij) and the static structure functions Sapik). Figs.(l) and (2) 
show these quantities in a 4.7% mixture, at a total density p = 0.02160A~'^, for the fLRiij) 
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correlation, in FHNC/0 approximation (i.e. we have neglected the elementary diagramaH^ 
The results for the SSF, with the fsR{r), differ mainly in the region of low-/c values, in 
agreement with the previous discussion. 

Table (1) shows the variational strenghts Z^p[k) of the one-phonon response for the 
same mixture and compares the results obtained with the SR and LR correlation functions 
at four momentum values, from q = 1.1 to 1.7A~^. The positions of the variational delta 
peaks, Uk, are also given. It has to be noticed that the Feynman spectrum overestimates 
the experimental data by at least lOK both in the maxon and roton regions. 

Table (2) provides the same quantities after the insertion of the 2PH states. Fig. (3) 
shows the '^He spectrum with the LRl correlation. The figure also compares the spectrum 
with pure ^He at po with the experimental results of ref.i (circles) in a X3 = 1.1% 
mixture at SVP and of ref.S (squares) for a X3 = 6.0% mixture. 

The changes in going from pure ^He to the mixture are clearly visible. This is mainly 
a density effect. In fact, we obtain similar results if the LR-ACA correlation is used. 
CBF perturbative corrections appear to be large and bring the maxon energy close to the 
experiments. The roton is not well described, as it is too shallow respect to the data. This 
feature is also present in the ^He case. As for the ^He spectrum, we believe that most of 
the discrepancy in this part of the spectrum has to be ascribed to the use of CA in the 
calculation of the CBF matrix elements. Moreover, contributions from higher order CBF 
pertubative diagrams are known to be important to correctly reproduce the roton minimum 
in pure ^HeS. However, the CA results show a change in the sign of the shift from mixture 
to pure system at g ~ 1.8A~^, in good agreement with the measured experimental value! 
at constant pressure {q ~ 1.9A~^). The boundaries of the Ip — Ih DSFs are related to the 
energies of the Ip — Ih state ep — eh- €p(^h) has been computed by the procedure of ref.0. 
However, because of the low ^He density, it turns out to be extremely close to the free Fermi 
gas spectrum 
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so, ef has been used in all the calculations. Perturbative corrections to ek may be 
computed in CBF. In the case of the ^He impurity, CBFPT provides an accurate evaluation 
of its spectrum0'0 if the decay of the impurity excited state (given by a correlated plane 
wave) into correlated IPH and 2PH states is considered. IPH states account for ~ 2/3 
of the difference between the experimental effective mass and the bare one, whereas 2PH 
states give the remainder. Because of the low density of the ^He component, it is reasonable 
to expect a similar behavior in the finite concentration mixture. It implies that we should 
insert in the CBFPT expansion the coupling between Ip — lh states and IPH and 2PH ones. 

Work along this line is in progress. Here we have used for e^^^ the CBFPT spectrum 
of the single impurity, obtained by extending to finite momenta the approach of ref.0 for 
the effective mass. The involved matrix elements have been computed in CA for the three- 
body distribution functions. CA gives m^{CA) ~ 1. 87/13 for the impurity, whereas the more 
realistic Superposition Approximation (SA) gives m^{SA) ~ 2.2m3. However, the SA k ^ 
matrix elements are much more involved than their CA counterparts, and their evaluation, 
together with a description of the method, will be presented subsequently^. Here, the effect 
of the missing effective mass has been estimated by simply scaling the CA spectrum as 
6^^^ = [m*iCA)/m*{exptM^^iCA). 

Fig. (4) shows the ^He spectra in different approximations, and compares them with the 
experimental data (circles from ref.i and squares from ref.S) and with the LP and LPM 
parametrizations given in the introduction, with parameters = 2.3m3 and 7 = 0.132A^. 
Even from the CA calculation, a deviation from the LP behavior clearly appears . The 
estimated CBF value of 7 in CA turns out to be -f{CA) = 0.052^^. We stress once more 
that we expect SA to provide a better description of the ^He spectrum behavior, as it 
correctly takes into account the core property of the system, requiring that the three-particle 
distribution functions vanish when any interparticle distance is lower than the radius of the 
repulsive core of the potential. 

Fig. (5) gives the Ip—lh DSF 5'^^^^^(g, u) at two momentum values {q = 1.3 and 1.7A~^), 
with the LRl correlation and using the CBF ^He spectrum. The two 33 DSF are very close 
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and dominant, becoming indistinguishable at higher momenta; the 44 component is always 
very small (notice that it has been amplified by a factor of 10 in the figure); the 34 part 
is negative and an order of magnitude larger than Sll'^^ in absolute value, contributing 
to decrease the total response mainly at low momenta. The free Fermi Gas DSF would 
be located to a larger energy with a lower peak strength, compatible with the fact that the 
correlated system has a ^He effective mass 2.3 times larger than the bare mass (at q = 1.3A~^ 
the FG peak position is c<j = 13. 67^ and the strength is Sl^'^^^FG) = .llSi^^^). In addition, 
as for the phonon DSF, the use of the SR and LR correlations does not alter appreciably 
the results shown in the figure. 

Table (3) shows the CBF values of the m„ Q,^(g) sum rules of the DSF's, for n = 0, 1, 
defined as: 

mn,af3iq) = J du uj"- S^p^q^uj)- (55) 
The exact DSF's satisfy the /—sum rules 

mi^apiq) = (56) 
and for n = one has mo,Q,/3(g) = Sa/siq)- 

The table gives also the variational values of the SSF's of eq.^ and the /—sum rules, 
mo and the SSF's, as well as mi and the /—sum rules are in good agreement for all g's 
in the 44 case. For 33, rriQ is reasonably close to the variational SSF, whereas the /—sum 
rule is underestimated (we recall that the SSF's have been computed with the variational 
ground states, not the exact one). The 34 m„ are not satisfying, especially for mi. In 
order to trace down the source of the differences we also show the contributions to the sums 
from the large-energy multiphonon tails of the DSF's, mi^Q{mPH). The tail provides a 
large part of the sum rules and, for ajS = 34, it is dominant. The effect is more clearly 
visible in mi, where the tails of the DSF's are multiplied by a large energy factor. As 
a consequence, a quantitatively correct estimate of mi^a/s would porbably require a more 
accurate evaluation of the multiphonon contributions as well as the insertion of higher order 
perturbative diagrams. 
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To evaluate the total scattering function S{q,uj), the DSF must be multiplied by the 
elementary cross sections and the concentrations of the species. In fig.s(6a,b) we give the 
partial CBF scattering functions (PSF): 

^44(g, uj) = — , , — -rSiAiq, w), (57) 

•^s^l?, ^) = ^ , w^r, ^^^ '^34(9, u;) (58) 
X4: + X3[a^ + CT3J 

and 

S33{q,u}) = — ——-r[cx3S33{q,uj) + alSl2{q,uj)]. (59) 

The LRl correlation has been used. The position and the strength of the phonon contri- 
bution to the PSF's are explicitely given. In the Ip— Ih sector, at the lower momentum, the 
33 PSF is strongly reduced by the 34 PSF, which practically disappears at g = 1.7A~^. The 
44 PSF is always negligible in this sector. In the phonon sector, the 44 PSF is the dominant 
one. The 33 component always results to be very small. The 34 PSF at g = 1.3A~^ slightly 
reduces the scattering function, while at g = 1.7A~^ increases it. 

In order to compare with the experimental scattering function, the theoretical PSF's 
have to be convoluted with the experimental broadening functions. As at these momentum 
transfers the phonon peak is still delta shaped, because there is no overlap with the multi- 
phonon background, we assume, in accordance with the authors of ref.i, that the width in 
energy of the low temperature results in that reference is entireky due to the instrumental 
resolution. For this reason we have convoluted the phonon peak with a gaussian having an 
average half maximum width of 1.3K. A gaussian with a width of 1.2id^ has been used for 
the Ip — Ih response. 

The convoluted total scattering functions are compared with the experimental results of 
ref.i in fig.s(7a,b,c) at g = 1.1, 1.5 and 1.7A~^, respectively, for the 4.7% mixture we have 
considered so far. 

At g = l.lA~^, both the position and the strength of the phonon branch are well de- 
scribed by our calculation. When approaching the roton minimun region, the agreement 
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worsens and we overestimate the experimental data. As discussed previously, we expect 
that the use of SA will improve the CBF description. 

An analogous analysis can be performed fort the Ip — Ih sector. The use of the CBF- 
CA spectrum slightly misses the location of the bump, well described in turn by a LPM 
parametrization, which is essentially a fit to the experimental data. We recall that the 
relevant difference between the LPM and the CBF-CA energies lies in the 7-parameter 
value, smaller by a factor ~ 0.4 in the latter case. A simple, quadratic LP parametrization 
with ml = 2.3ms seems to be ruled out. 

The ^He scattering function 5*3(^,0;), defined as 

S3fe-)-^^^|±^^%-). (60) 

and the function 5*3(^,0;), given by 

in the Ip — Ih sector, are given in fig. (8). The figures contains also the convolution of 
Ss{q,uj) with the experimental broadening function and the experimental results of ref.i, 
at g = 1.3A~^. S3 and ^3 are identical if S^^~^'^ = Sl^~^'^ = 0. So, their differences are 
basically a measure of the importance of the 34 contribution (the 44 one being negligible). 
Our results show a large suppression of 5*33^^'^ due to 5*34"^'^, which brings the CBF response 
much closer to the experiments. 

5. CONCLUSIONS 

Correlated Basis Perturbation Theory has been used to microscopically compute the 
scattering function in a X3 = 4.7% ^He-^He mixture at T=0. The theory has allowed for 
explicitely separating the different contributions to the response and for semiquantitatively 
assessing the relevance of the 34 component. In the Ip — Ih region, the 6*33 response is 
sizeably reduced by 5*34 up to g ~ 1.5A~\ whereas 5*44 is always negligible. A similar effect, 
even if smaller in magnitude, is present in the phonon-roton sector, where the dominant 6*44 
is only sligthly modified by 6*34. 



The responses have been computed by inserting correlated Ip — Ih and 1- and 2-phonon 
intermediate states. Also the possible decay of 1-phonon into 2-phonon states has been 
estimated in boson-boson approximation and using the Convolution Approximation for the 
three-body distribution functions. 

The microscopic quasiparticle ^He spectrum clearly shows a deviation from the simple 
LP form. The spectrum has been actually computed for the single impurity problem, but 
we do not believe that its evaluation in the low concentration mixture will dramatically 
change our findings. In particular, a deviation from LP was advocated in ref.i to explain 
the experimental Ip — Ih response, in contrast with a possible large change of the ^He 
effective mass in mixture (from mg = 2.3m3 at X3 = to mg = 2.9m3 at X3 = 4.7%). The 
CBF spectrum still does not reproduce fully quantitatively the data, and a more accurate 
calculation is needed. 

The '^He excitation spectra in the phonon-roton branch of the pure system and the 
mixture at SVP have been compared. The shift between the two excitations appears to 
be due to the change in density. CBF gives a good description of the maxon region, but 
overestimates the roton, even if it gives an almost correct g-value for the change of sign of 
the shift . 

The CBF scattering function at low momenta gives a reasonable description of the scat- 
tering data (both for the position and strength). The agreement worsens as q increases. The 
peaks are located at too a large energy and their strength is overestimated. We believe that 
the reason of this lies in the approximations made to compute the decay of IPH states into 
2PH and in the lack of higher intermediate states, which become more and more important 
as the momentum increases. In particular, the Ip—lh sector does not include two probably 
relevant contributions: the decays of Ip — Ih states into (1) 2p — 2h and (2) IPH states. 
The former adds large energy tails to the Ip — lh bump reducing its strength, and the latter 
is known to be responsible for a large part of the ^He effective mass. Our CBF calculation 
includes the real part of the Ip — Ih into IPH decay but does not consider its imaginary 
part. 
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The importance of the 34 contribution to the total scattering function is especially visible 
in the ^He scattering function in the Ip — Ih region, where its introduction reduces the 
response by a factor ~ 0.6 at g = 1.3A^^. 

More work is clearly needed in order to give a fully quantitative description of both the 
excitations and the responses of the Helium mixtures. However, from our results, we believe 
that CBF is a promising theory in view of achieving this goal. 
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FIGURES 

FIG. 1. Radial distribution functions for the mixture (see text). The sohd hne gives 544, the 
dashed hne 5134 and the dotted one is 5133. 

FIG. 2. Static structure functions for the mixture (see text). The solid line gives 5*44, the 
dashed line S34 and the dotted one is 833. 

FIG. 3. ^He excitation spectrum in the mixture (crosses) and in the pure system (dashed 
line). The upper curves are the Feynam spectra. Squares and circles are mixture experimental 
data (See text). 

FIG. 4. ^He excitation spectrum. The dashed line gives e^^^ (CA), the solid line is e^^^^', 
also shown are the LP, LPM and free (F) spectra. Squares and circles are the experimental data 
(See text). 

FIG. 5. Ip — Ih DSF at q = 1.3 and 1.7 A~^. The continuous line gives S'33,the dot-dashed 
5*33, the dotted 10 x 6*44 and the dashed 5*34. 

FIG. 6. CBF Partial scattering functions at q = 1.3 (6a) and 1.7A~^ (6b). Continuous line 
S33, dashed line S^^, dotted S'44. The PH-q;/3 numbers are the strengths of the phonon PSF, 
located at LOq. 

FIG. 7. Total scattering functions at g = 1.1 (7a), q = 1.5 (7b) and 1.7A~^ (7c) (sohd lines). 
Also shown are the Ip — Ih responses with the LP (dotted lines) and LPM (dashed lines) spectra 
and the experimental data (crosses). 

FIG. 8. Ip — Ih ^He scattering function at q = 1.3A~^. The solid line is S3{q,u!), the dashed 
line is S3{q,Lo), the dotted line is the experimental convolution of 5*3(5,0;) and crosses are the 
experimental data. 
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TABLES 

TABLE I. Variational strenghts and positions of the one phonon DSF responses with different 

correlations (see text), q in A^^ and ujq in K 





Q 


LOg 




^34 


yv 
^33 


SR 


1.1 


20.44 


0.356 


-0.146 


0.060 


LR 




20.77 


0.361 


-0.145 


0.058 


LRl 




20.75 


0.362 


-0.147 


0.060 


SR 


1.3 


20.74 


0.491 


-0.115 


0.027 


LR 




21.68 


0.461 


-0.122 


0.032 


LRl 




21.71 


0.460 


-0.122 


0.032 


SR 


1.5 


19.91 


0.681 


-0.072 


0.008 


LR 




19.99 


0.686 


-0.071 


0.007 


LRl 




20.02 


0.685 


-0.069 


0.007 


SR 


1.7 


19.20 


0.908 


-0.021 


0.000 


LR 




19.24 


0.899 


-0.023 


0.001 


LRl 




19.26 


0.899 


-0.019 


0.000 
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TABLE II. CBF strenghts and positions of the one phonon DSF responses with different cor- 



relations (see 


text), q in A ^ 


and LOq in K 










q 




Z44 




Z33 


SR 


1.1 


13.73 


0.275 


-0.066 


0.016 


LR 




13.69 


0.272 


-0.068 


0.017 


LRl 




13.66 


0.272 


-0.068 


0.017 


SR 


1.3 


14.01 


0.390 


-0.045 


0.005 


LR 




14.27 


0.367 


-0.047 


0.006 


LRl 




14.25 


0.366 


-0.044 


0.005 


SR 


1.5 


13.84 


0.559 


-0.014 


0.000 


LR 




13.79 


0.558 


-0.014 


0.000 


LRl 




13.83 


0.557 


-0.011 


0.000 


SR 


1.7 


13.94 


0.766 


0.024 


0.001 


LR 




13.94 


0.766 


0.024 


0.001 


LRl 




13.96 


0.765 


0.027 


0.001 
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TABLE III. CBF sum rules mo and mi with the LRl correlation and CBF ^He spectrum, 
variational SSF's and /-sum rules, q in and energies in K. In parentheses are shown the 
multiphonon tail contributions 



wo(rnPn) 



miimPII) 



f-SR 



44 
34 
33 

44 
34 
33 

44 
34 
33 

44 
34 
33 



1.1 



1.3 



1.5 



1.7 



0.37(0.08) 
-0.29(-0.08) 
1.10(0.13) 

0.49(0.10) 
-0.24(-0.08) 
1.05(0.09) 

0.68(0.12) 
-0.14(-0.06) 
1.00(0.05) 

0.91(0.14) 
-0.04(-0.04) 
0.96(0.03) 



7.38(3.58) 
-4.80(-3.28) 
8.88(4.67) 

10.27(4.97) 
-4.53(-3.22) 
8.97(3.47) 

13.59(5.86) 
-3.18(-2.53) 
9.15(2.14) 

17.40(6.72) 
-1.66(-1.86) 
9.93(1.25) 



0.36 
-0.15 
0.98 

0.46 
-0.12 
0.98 

0.68 
-0.07 
0.99 

0.90 
-0.02 
1.00 



7.29 
0.0 
9.73 

10.19 
0.0 
13.59 

13.56 
0.0 
18.09 

17.42 
0.0 
23.24 
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